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Abstract 

The pseudo-perturbative shifted - I expansion technique PSLET is shown 
applicable in the non-Hermitian VT - symmetric context. The construction of 
bound states for several VT - symmetric potentials is presented, with special 
attention paid to V(r) = ir 3 — ayir oscillators. 



I. INTRODUCTION 



In their recent studies Dorey, Dunning and Tateo (DDT) [1] have considered the mani- 
festly non-Hermitian Schrodinger equation, in h = 2m = 1 units, 



d 2 £{£+!) r- . „ 

-r^7 H t, aVir + ir° 

dr A r z 



ipk,e( r ) = E k/ ip k/ (x). (1) 



They have rigorously proved that the spectrum E kj g is real and discrete in the domain of 
the sufficiently large angular momenta, 



£ > max 



i(2«-7),-i 



(2) 



This inspired our subsequent study of this model [2] where we have shown that in the strong 
coupling regime with £ >> 1, the low lying DDT bound states may be very well approximated 
by the harmonic oscillators. At the same time we have noticed that the quality of such 
an asymptotic approximation may deteriorate quite significantly with both the increase of 
excitation k and/or with the decrease of £. 

Such a situation is, obviously, challenging. Firstly, our study [2] revealed that the mani- 
fest non-Hermit icity of the models of the type (1) leads to the reliable leading order approx- 
imation only after we select our harmonic oscillator approximant as lying very far from the 
real axis (i.e., from the Hermitian regime). Such a recipe is, apparently, deeply incompatible 
with a smooth modification of the traditional zero-order approximants occurring in current 
Hermitian \ j£ recipes (cf. a small sample of some references in [3]). At the same time, the 
smallness of l/£ still supports the feeling that the similar perturbation techniques should 
prove efficient after their appropriate modification. 

This observation offered a sufficiently strong motivation for our continued interest in the 
complex, non-Hermitian model (1) which may be understood as a characteristic represen- 
tative of a very broad class of the so called pseudo-Hermitian models with real spectra, the 
analyses of which became very popular in the recent literature [4] - [18]. Within this class, 
the strong-coupling version of DDT oscillators (1) with a ^> 1 forms a particularly suitable 
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testing ground as it combines the necessary reality of its spectrum with the smallness of 
the inverse quantity Moreover, the phenomenologically appealing non-Hermitian mod- 
els like (1) are rarely solvable in closed form so that the presence of a "universal" small 
parameter l/£ <C 1 offers one of not too many ways towards their systematic approximate 
solution. 

In our present paper II we intend to discuss such a possibility in more detail. 

II. FRAMEWORK 

The first stages of interest in the non-Hermitian oscillators (1) date back to an old paper 
by Caliceti et al [5]. It studied the imaginary cubic problem in the context of perturbation 
theory and, more than twenty years ago, it offered the first rigorous explanation why the 
spectrum in such a model may be real and discrete. In the literature, this result has been 
quoted as a mathematical curiosity [6] and only many years later, its possible relevance 
in physics re-emerged and has been emphasized [7]. This initiated an extensive discussion 
which resulted in the proposal of the so called VT symmetric quantum mechanics by Bender 
et al [8]. 

The key idea of the new formalism lies in the empirical observation that the (phenomeno- 
logically desirable) existence of the real spectrum need not necessarily be attributed to the 
Hermiticity of the Hamiltonian. The current Hermiticity assumption H = if' is replaced 
by the mere VT symmetry H = = VTHVT . Here, V denotes the parity (VxV = —x) 
while the anti-linear operator T mimics the time reflection [TiT = —i). It is easy to verify 
that example (1) exhibits such a type of symmetry [2] and may serve as an elementary 
illustration of the latter extension of quantum mechanics. 

The Bender's and Boettcher's conjecture that H = H$ may imply ImE = is fragile. 
The extent as well as limitations of its validity are most easily analyzed in the language of 
linear algebra using the biorthogonal bases [4,9] and/or exactly solvable Hamiltonians [10]. 
Nevertheless, the relevance of many unsolvable oscillators originates from their applicability 
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in physics [11] and field theory [12]. In such a setting it is necessary to develop and test also 
some efficient approximation methods. New and intensive studies employed the ideas of the 
strong-coupling expansions [13] as well as the complex version of WKB [14], Hill determi- 
nants and Fourier transformation [15], functional analysis [16], variational and truncation 
techniques [17], and linear programming [18]. 

In what follows we intend to use the method based on the smallness of the inverse an- 
gular momentum parameter l/£. Various versions [19] of such an approach are available 
for Hermitian models where the combination of the central repulsive core £(£ + l)/r 2 with 
a confining (i.e., asymptotically growing) interaction V(r) forms a practical effective po- 
tential V p {r) which possesses a pronounced minimum. Near such a minimum the shape of 
the potential is naturally fitted by the elementary and solvable harmonic oscillator well. 
Corrections can be then evaluated in an unambiguous and systematic manner [20]. 

As long as we intend to move to the complex plane, the leading-order approximation 
may become non-unique. One finds several different complex and/or real minima of V p (r) 
even in our oversimplified examples (1) [2]. For all the similar non-Hermitian Hamiltonians, 
even the most sophisticated forms of the perturbation expansions in the powers of our small 
parameter l/£ lose their intuitive background and deserve careful new tests, therefore. 

Some of the related reopened questions will be clarified by our forthcoming considerations 
inspired, basically, by the pseudo-perturbative shifted— £ expansion technique (PSLET) in 
its form designed for the standard Hermitian Hamiltonians and described, say, by Mustafa 
and Odeh in ref. [20]. In this recipe, the shifts of £ were admitted as suitable optional 
auxiliary parameters while the use of the prefix "pseudo-" just indicated that £ itself is an 
artificial, kinematical parameter rather than a genuine dynamical coupling. 

III. VT -SYMMETRIC PSLET RECIPE 

As we already mentioned, one of the first VT — symmetric models with an immediate 
impact on physics has been the Buslaev's and Grecchi's quartic anharmonic oscillator [7] 



4 



described by the radial Schrodinger equation in d— dimensional space, 

XkA r ) = E k,e XkA r )- ( 3 ) 

In this model they shifted the coordinate axis to the complex plane, r = £ — ic with a 
constant S(r) = — c < and variable 5ft(r) = £ G (—00,00). They also required that 
Xk,e( r ) e -^2(^00, 00) at all the partial waves id = £ + (d — 3)/2 and dimensions d > 2. 

This example may find the various sophisticated generalizations some of which will be 
also mentioned in due course in what follows. For example, a £— dependent shift c = c(£) 
may be needed both for the exactly solvable Coulombic model of ref.[21] and for all the 
more general and purely numerically tractable potentials V(r) ~ — (ir) N of ref.[8] with the 
positive exponents N > 3. Fortunately, the transition to c = c(£) remains particularly 
elementary, being mediated by the mere change of the variable in eq.(3) within this class 
(c.f., e.g., ref. [21] for an explicit illustration). Another remark might concern the assumption 
that the wave functions are square integrable. For the most elementary VT— symmetric 
Hamiltonians this assumption seems very natural but in certain more sophisticated models 
its use may require a more careful analysis as presented, e.g. in refs. [4] and [16]. 

The practical experience with the Hermitian version of the pseudo-perturbation shifted - 
£ expansion technique of Mustafa and Odeh [20] may serve as a key inspiration of an appro- 
priate complexified new PSLET recipe. Firstly we notice a formal equivalence between the 
assumed smallness of our parameter 1 / id ~ and of its ( arbitrarily) shifted form. Thus 
we introduce a new symbol I — td~ (3 and, simultaneously, move and rescale our coordinates 
r — > x to 



dr 2 



+ ''"', + 1 > +y( r) 



x = P /2 (r -r ). (4) 

Here r is an arbitrary point, with its particular value to be determined later. Equation (3) 
thus becomes 

T d 2 I 2 + {2(3 + 1)1 + (3 {(3 + 1) I 2 T „ . ..1 T . . / / \ \ , c x 

I ^ + r g Mro \- 1/2)+ \ ]2 + Q "(*(r))} = E k ^ kA x{r)l (5) 



where Q is a constant that scales the potential V at large-£d limit and is set, for any specific 
choice of id and k, equal to P at the end of calculation. Expansions about this point, x = 
(i.e. r — ro), yield 



[x/(r„P/2) + l] 2 „.„ 



I V ™) = E (^) ^ '- ( "- 4V2 < 7 » 
It is also convenient to expand E^/ as 

oo 

E k ,e = E *V ( g ) 

n=-2 

Of course one may consider also energy coefficient of half-entire power of / in (8) but all 

these coefficients vanish ( c.f., e.g., ref. [19,20]). Equation (5), therefore, reads 



dX n=0 n=0 



+/3 (/? + 1) E T « ^ r " (n+2)/2 >*m(*) 

n=0 

oo ^ 

E 4? r_(n+1) *MW (9) 



,n=-2 



where 



Equation (9) is to be compared with the non-Hermitian VT - symmetrized perturbed har- 
monic oscillator in the one dimensional Schrodinger equation 



— 2 + — [y - %c ) + £ o + "u/ - ic) 



$ fe (y) = A fe $ fe (y), 



(11) 



cfr/ 2 ' 4 

where P(y — ic) is a complexified perturbation - like term and eo is obviously a constant. 
Such a comparison implies 
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e = B I + (2(3 + 1)T + (3((3 + 1)T / 1 



n=0 



= B l + [(2P + 1)T + (2k + l)^]+i 0(0 + 1)T + A< 0) I + £ Ajf _1) / 



2 J / 



ra=2 



(12) 



n=l 



The first two dominant terms are obvious 



F (- 2 ) _ 1 , V(r ) 



(13) 



r 5 



and with appropriate rearrangements we obtain 



< = M ^ + A«°», (15) 







E$ = A<»> ; n > 1. (16) 



Here r is chosen to minimize p , i. e. 



dr dr^ 
Equation (13) in turn gives, with (£ d — (3) 2 = Q, 



= and ,T > 0. (17) 



i^,.^). (18) 

Consequently, B 1 ( = I E kl ) adds a constant to the energy eigenvalues and B\ = 0. 
The next leading correction to the energy series, I E^\ consists of a constant term 
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and the exact eigenvalues of the unperturbed one-dimensional harmonic oscillator poten- 
tial uj 2 x 2 /4 (= B 2 x 2 ), where 



0<u, = u>W = ± ^Q; Q = j3+ r ° V " {ro) 



Evidently, equation (19) implies that r$ is either pure real, uo = uj( + \ or pure imaginary, 
uj = uj(~\ Next, the shifting parameter (3 is determined by choosing T E^~p = 0. That is 



(19) 



/? = /? (±) = ~[i±(2 fc + i) n], 



(20) 



where (3 = ft +s> for r pure real and (3 = ft ) for r pure imaginary. Then equation (9) 
reduces to 



dx 2 



+ J2 v(n) J ~ n 



/2 



n=0 



n=-l 



(21) 



with 



,(°\x) = B 2 x 2 + (2[3+l) T , 



(22) 



and for n > 1 



, {n \x) = B n+2 x n+2 + (2(3 + l)T n x n + (3((3 + l) T n _ 2 x n ~ 



Equation (21) upon setting the wave functions 



^k,e( x ) = F k,e( x ) exp(U k/ (x)) 



(23) 



readily transforms into the Riccati-type equation: 
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F k A x ) 



CO 

- (u'U*) + u' k ,i(x)u' k M) + I> (n) 0#- 



■n/2 



n=0 



n=0 



-2F' k /x)U' k /x)-F' k [ e (x) = 0, 



where the primes denote derivatives with respect to x. It is evident that this equation admits 
solutions (c.f., e.g., ref.[20]) of the form 



?1=0 



n=0 



with 



00 k— 1 



n=0 p=0 



Ui n \x) 



n+l 
m=0 



2m-l 



! D(),n,k — 0, 



n+l 

^k \ X ) — ^ ^ m '™: fc X2m - 

m=0 

Obviously, equating the coefficients of the same powers of I and x ( for each k), respectively, 
one can calculate the energy eigenvalues and eigenf unctions ( following the uniqueness of 
power series representation) from the knowledge of C mi n : k, D mjn: k, and a^J. in a hierarchical 
manner. 

In order to test the performance of our strategy, let us first apply it to the two trivial, 
exactly solvable VT- symmetric examples. 
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IV. AN ELEMENTARY ILLUSTRATION OF THE RECIPE 



A. VT— symmetric Coulomb 

Using the potential V(r) = iA/r (where A is a real coupling constant) in the above 
"PT— symmetric PSLET setting, one reveals the leading-order energy approximation 

= -2 + -■ (24) 

The unique minimum at r = 2iP/A occurs in the upper-half of the complex plane. In this 
case Q — 1, (3 = (3^ = k, the leading energy term reads 



p (_ 2) _ iA _ A 2 



and higher-order corrections vanish identically. Therefore, the total energy is 

where n = k + £ + lis the principle quantum number. Evidently, the degeneracy associated 
with ordinary ( Hermitian) Coulomb energies E n = —A 2 /(2n) 2 is now lifted upon the com- 
plexification of, say, the dielectric constant embedded in A. Moreover, the phenomenon of 
flown away states at k = id emerge, of course if they do exist at all ( i. e. the probability 
of finding such states is presumably zero, the proof of which is already beyond our current 
methodical proposal). Therefore, for each k— state there is an state to fly away. 

Next, let us replace the central-like repulsive/attractive core through the transformation 
£d(£d + 1) c? ~~ ^-/^i i- e - — — 1/2 + q\a \, with q — ±1 denoting quasi-parity, and recast 



(25) as 



A 2 

(2k + 1 - 2q\a \ 



Which is indeed the exact result obtained by Znojil and Levai [21]. Equation (27) implies 
that even- quasi-parity, q = +1, states with k = \a \ — 1/2 fly away and disappear from 
the spectrum. Nevertheless, quasi-parity- oscillations are now manifested by energy levels 
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crossings. That is, a state-A; with even-quasi-parity crosses with a state-A;' with odd-quasi- 
parity when \a Q \ — (k — k) /2. However, when a Q = the central-like core becomes attractive 
and the corresponding states cease to perform quasi-parity-oscilations. For more details on 
the result (27) the reader may refer to Znojil and Levai [21]. 

B. VT— symmetric harmonic oscillator V(r) = r 2 

For this potential the leading energy term reads 

'X^S + ro 2 , (28) 
'o 

and supports four eligible minima ( all satisfy our conditions in (17)) obtained through 
Tq = P as r = ± i\ I I 1 / 2 and r = ± \ I \ 1 / 2 . In this case Q = 2, 

/3 = = - (2k + 3/2) (29) 

for r = ± | I I 1 / 2 , and 

= /3H = (2fc + 1/2) (30) 
for r = ± i | f I 1 / 2 . Whilst the former (29) yields 

PE { - e 2) = 2rl = Ak + 2i d + 3, (31) 

the latter (30) yields 

pEj-p = 2rl = Ak + l- 2£ d . (32) 

In both cases (3 = (3^ the higher-order corrections vanish identically. Yet, one would 
combine (31) and (32) by the superscript (±) and cast 

£g } =4fc + 2±(24i + l). (33) 

Therefore, the VT— symmetric oscillator is exactly solvable, by our recipe, and its 
spectrum, non-equidistant in general, exhibits some unusual features ( cf. [22] for more 
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details). However, it should be noted that for the one-dimensional oscillator ( where ^ = — 1, 
and 0, even and odd parity, respectively) equation (33) implies (I) E^/2 = 2k + 1/2, 
£H/2 = 2k + 3/2, for £ d = -1, and (II) E^/2 = 2k + 3/2, £?H/2 = 2A; + 1/2, for £ d = 0. 
Which can be combined together to form the exact well known result 

E N = 2N+1; N = 0,1,2,.... (34) 

with a new, redefined quantum number N. 



V. APPLICATION: VT -SYMMETRIC DDT OSCILLATORS 



In our VT— symmetric Schrodinger equation (1) with the practical effective potential 



V p (r) = ^"[^ -a/^ + tr 3 (35) 



the general solutions themselves are analytic functions of r (c.f., e.g., [6]). We may construct 
them in the complex plane which is cut, say, from the origin upwards. This means that 
r = £ exp(i(p) with the length £ e (0, oo) and with the span of the angle (p e (— 37r/2, it/2). 
Compact accounts of the related mathematics may be found in Bender et al [8]. 

Let us proceed with our PT-PSLET and search for the minimum/minima of our leading 
energy term for the ZTDr-oscillators (1) 

72 

PE { k - 2) = — - a^fTFo + % r 3 . (36) 
r o 

Evidently, condition (17) yields 

I"1 _ 9 _1 

0. (37) 



r 5 + % 



1 2 

-a ro 2 \fTr~Q + - I 
6 3 



Obviously, a closed form solution for this equation is hard to find ( if it exists at all) and 
one has to appeal to numerical techniques to solve for r . 

A priori, it is convenient to do some elementary analyzes, in the vicinity of the extremes 
of a ( mandated by condition (2)), and distinguish between the two different domains of 
a. For this purpose let us denote (2P)/3 = G 2 , re-scale r = — i\G 2 ^\ p and abbreviate 
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a = 5VQ G 2 with < S < 1. This gives the following new algebraic transparent form of 
our implicit definition of the minimum/minima in (37), 

1 - Z = 5 ^1 ; . Z = p\ (38) 

In the weak-coupling domain, vanishing 5 ~ 0, equation (38) becomes trivial (1 — Z — 0). 
It is easy to verify that (36), with 5^0, has a unique absolute minimum at 



Z = 1 =^p = 1 =^r = - i\G 2/5 \ , 5 = 0. 



(39) 



In the strong-coupling regime, 5 ~ 1, equation (38) yields Z = 2/3. 

At this point, one may choose to work with (3 — ( i.e. -E^^ 7^ 0) an d obtain the 
leading (zeroth)-order approximation 

15" 



Pe { - 2) 



q6\ V 5 

~Z^ 



5Z- 



and, with 



U R 

T 2 



' z 



r = -iG 2/5 Z 1/5 , 



the first-order correction 



(-i) 

kl 



2 

G 
Z 2 



G 



' +(2* + l)£ 



L r 



1/5 



15 (1 + Z) 



(2k + 1) 



Consequently, the energy series (8) reads, up to the first-order correction 



E. 



Nevertheless, one may choose to work with (3 — (3^ ^ ( i.e. e\, p = 0) and obtain 



(40) 



(41) 



(42) 



(3 = (3 



(-) - — 



1 - ( 2^ + 1)^5 (l + Z)/2 



(43) 



Thus the zeroth-order approximation yields 
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In Figure I we plot the energies of (44) vs t G (—5, 5) at different values of Z = Z(5) G 
(2/3, 1). Obviously, our results show that even with £ < £ddt («) the spectrum remains real 
and discrete. Moreover, once we replace £(£ + 1) — > a 2 Q — 1/4, i.e. i — >■ —1/2 + q\a Q \ with 
q — ±1 denoting quasi-parity, quasi-parity-oscillations are manifested by the unavoidable 
energy levels crossings, see Figure II. 

Table I shows that our results from Eqs.(42) and (44) compare satisfactorily with those 
obtained by Znojil, Gemperle and Mustafa [2], via direct variable representation (DVR). We 
may mention that even in the domain of not too large £, the difference between the exact 
and approximate energies remain small, of order of « 0.05% from Eq.(42), with (3 = 0, 
and ~ 0.2% from Eq.(44), with (3 = (3^ ^ 0, for the ground state. Such a prediction 
should not mislead us in connection with the related convergence/divergence of our energy 
series (8), which is in fact the genuine test of our present VT— symmetric PSLET formulae. 
The energy series (8) with (3 = (3^ ^ converge more rapidly than it does with (3 — 0. 
Nevertheless, our leading energy term remains the benchmark for testing the reality and 
discreteness of the energy spectrum. 

In Table II we compare our results ( using (3 = (3^ ^ 0, hereinafter, numerically solve 
for r and following the procedure of section III) for (1) with a = using the first ten-terms 
of (8) and the corresponding Pade approximant, again with those from the DVR approach. 
They are in almost exact accord. Hereby, we may emphasize that the digital-precision 
enhances for larger I ( smaller l/l ) values, where the energy series (8) and the related Pade 
approximants stabilize more rapidly. 

Extending the recipe of our test beyond the weak-coupling regime 5 ~ (i.e. a w 0) 
we show, in table III, the energy dependence on the no n- vanishing a. Evidently, the digital- 
precision of our VT - PSLET recipe reappears to be /-dependent and almost a - independent. 
In table IV we witness that the leading energy approximation inherits a substantial amount 
of the total energy documenting, on the computational and practical methodical side, the 
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usefulness of our pseudo-perturbation recipe beyond its promise of being quite handy. Yet, 
a broad range of a is considered including the domain of negative values, safely protected 
against any possible spontaneous VT -symmetry breaking. 

VI. SUMMARY 

Our purpose was to find a suitable perturbation series like expansion of the bound states. 
We have started from the observation that the physical consistency of the model (1) (i.e., 
the reality of its spectrum) is characterized by the presence of a strongly repulsive/attractive 
core in the potential V p (r). This is slightly counterintuitive since the phenomenologically 
useful values of £ are usually small and only the first few lowest angular momenta are relevant 
in the Hermitian Schrodinger equations with central symmetry. 

Only exceptionally, very high partial waves are really needed for phenomenological pur- 
poses (say, in nuclear physics [23]). The strong repulsion is required there, first of all, due 
to its significant phenomenological relevance and in spite of the formal difficulties. 

Fortunately, efficient £ ^> 1 approximation techniques already exist for the latter partic- 
ular realistic Hermitian models. They have been developed by many authors (cf., e.g., their 
concise review in [24]). Their thorough and critical tests are amply available but similar 
studies were still missing in the non-Hermitian context. 

Our present purpose was to fill the gap at least partially. We have paid thorough attention 
to the first few open problems related, e.g., to the possible complex deformation of the axis of 
coordinates. Our thorough study of a few particular VT symmetric examples revealed that 
the transition to the non-Hermitian models is unexpectedly smooth. We did not encounter 
any serious difficulties, in spite of many apparent obstacles as mentioned in Section II (e.g., 
an enormous ambiguity of the choice of the most suitable zero-order approximation). 

In this way, our present study confirmed that the angular momentum (or dimension) 
parameter t in the "strongly spiked" domain where \l\ >> 1 offers its formal re- interpretation 
and introduction of an artificial perturbation- like parameter 1/7 which may serve as a guide 
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to the interpretation of many effective potentials as a suitably chosen solvable (harmonic 
oscillator) zero-order approximation followed by the systematically constructed corrections 
which prove obtainable quite easily. 
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FIGURE CAPTIONS 

Figure 1. Cubic oscillator (35) eigenenergies (— E) in (44) vs I for k = and | i \ < 5 
at different values of Z = Z(5) G (2/3, 1). 

Figure 2. Cubic oscillator (35) eigenenergies (— E) in (44) vs \a Q \ for k = 0, I = 
—1/2 + q\a \, and different values of Z = Z{8) e (2/3, 1) at even- and odd-quasi-parities. 
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TABLES 

TABLE I. Comparison of the energy levels for model (1) ( with a = 0). The benchmark, 
numerically exact DVR values are cited from Znojil, Gemperle and Mustafa [2]. 



t 


k 


DVR 


Eq.(42) 


Eq.(44) 


5 





-11.52191 


-11.517 


-11.542 




1 


-4.56482 


-4.260 


-4.900 




2 


1.87017 


2.997 


-0.109 


10 





-28.76552 


-28.762 


-28.776 




1 


-20.59867 


-20.426 


-20.756 




2 


-12.70640 


-12.090 


-13.230 




3 


-5.11663 


-3.754 


-6.380 




4 


2.14032 


4.58 


-0.727 


20 





-68.72646 


-68.724 


-68.733 




1 


-59.24706 


-59.149 


-59.330 




2 


-49.91773 


-49.574 


-50.171 




3 


-40.74589 


-39.998 


-41.283 




4 


-31.73951 


-30.423 


-32.705 




5 


-22.90712 


-20.847 


-24.489 




6 


-14.25769 


-11.272 


-16.713 




7 


-5.80054 


-1.696 


-9.512 




8 


2.45491 


7.879 


-3.172 


50 





-211.13555 


-211.134 


-211.138 




1 


-199.68009 


-199.633 


-199.718 




2 


-188.29459 


-188.132 


-188.406 




3 


-176.98040 


-176.631 


-177.206 




4 


-165.73889 


-165.129 


-166.123 




5 


-154.57149 


-153.628 


-155.161 
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6 -143.47967 -142.127 -144.328 

7 -132.46494 -130.626 -133.628 

8 -121.52886 -119.124 -123.069 
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TABLE II. Same as table I with PT-PSLET results from the first ten-terms of (8) and the 
corresponding Pade approximant. 



k 


I 


DVR 


PT-PSLET 


Pade 


r» 





5 


-11.52191 


-11.52191 


-11.52191336 


4.4 




10 


-28.76552 


-28.76552 


-28.76552178 


9.4 




20 


-68.72646 


-68.72646 


-68.72645928 


19.4 




50 


-211.13555 


-211.135548 


-211.13554785 


49.4 


1 


5 


-4.56482 


-4.565 


-4.564813 


2.2 




10 


-20.59867 


-20.598669 


-20.59866911 


7.2 




20 


-59.24706 


-59.24705556 


-59.247055572 


17.2 




50 


-199.68009 


-199.68008657 


-199.6800865747 


47.2 


2 


10 


-12.70640 


-12.7065 


-12.7964017 


4.9 




20 


-49.91773 


-49.917727 


-49.917727248 


14.9 




50 


-188.29459 


-188.294591 


-188.29459127507 


44.9 


3 


10 


-5.11663 


-9.388 


-5.1166 


2.7 




20 


-40.74589 


-40.74589 


-40.74589026 


12.7 




50 


-176.98040 


-176.98039968 


-176.98039968 


42.7 
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TABLE III. Comparison of the energies for model (1) with I = 10 and different values of a. 



a 
20 



DVR 



PT-PSLET 



Pade 




1 
2 



-58.62190 
-49.83626 
-41.29014 



-58.6219026 
-49.836258 
-41.2905 



-58.621902667 
-49.8362579 
-41.29014 



10 




1 

2 



-43.85223 
-35.40717 
-27.23003 



-43.85223324 
-35.407174 
-27.23010 



-43.852233235 
-35.40717408 
-27.23003 




1 

2 



-28.76552 
-20.59867 
-12.70640 



-28.76552 
-20.59867 
-12.70653 



-28.76551777 
-20.598669108 
-12.7064017 



-10 




1 
2 



-13.35529 
-5.40717 
2.27617 



-13.3552878 
-5.40717 
2.08 



-13.355287796 
-5.4071736 
2.27617 



-20 




1 
2 



2.38131 

10.16462 

17.70339 



2.381306 
10.1646 



2.38130557 
10.1646191 
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TABLE IV. Energies for model (1) with I = and different values of a, d, k. 



dak Leading term PT-PSLET Pade I « 

3 -5 3.07 2.903 2.881 1 

1 -1.65 -1.58 -1.579391 3.2 

2 -7.962 -7.652 -7.65505 5.4 
-10 8.23 7.82 7.820 1.3 

1 3.98 3.803 3.80992 3.5 

2 -1.87 -1.80 -1.798903 5.7 
-15 13.95 13.38 13.38419 1.6 

1 9.95 9.52 9.542833 3.8 

2 4.45 4.27 4.273163 6.0 
1 10 -12.46 -12.471 -12.471 1.4 

1 -20.65 -20.304 -20.3040 3.5 

2 -29.00 028.341 -28.34038 5.7 
5 -8.12 -8.082 -8.082 1.5 

1 -15.38 -15.058 -15.05795 3.6 

2 -23.15 -22.569 -22.56820 5.9 
-5 1.57 1.5393 1.5393 1.8 

1 -4.16 -4.055 -4.055465 4.0 

2 -10.94 -10.63 -10.6341055 6.3 
-15 12.99 12.754 12.75391 2.2 

1 8.06 7.84 7.84006 4.6 

2 1.99 1.93 1.9269341 6.8 
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Figure 1 




Figure 2 
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